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LAMINAR  BOUNDARY  LAYER  OVER  FLAT  PLATE  IN  A  FLOW  HAVING 
CIRCULAR  STREAMLINES 
By  Artur  Mager  and  Arthur  G.  Hansen 


SUMMARY 

The  laminar-houndary-layer  development  over  a  semi-infinite  flat 
plate  placed  in  a  flow  with  concentric  circular  streamlines  was  investi¬ 
gated  with  the  limitation  of  small  total  turning  of  the  main- stream 
flow.  The  shape  of  the  velocity  profiles  in  the  direction  of  the  main¬ 
stream  flow  and  perpendicular  to  it  was  analytically  determined  for  an 
incompressible  flow  and  a  compressible  flow  with  Prandtl  number  equal 
to  1. 


,  The  boundary-layer  thickness  was  shown  to  be  proportional  to  the 

square  root  of  the  distance  from  the  leading  edge  of  the  plate  when 
measured  along  the  streamline  of  the  main-stream  flow.  The  deflection 
*  of  the  boundary-layer  flow  at  the  plate  surface  from  the  direction  of 

a  circular  streamline  in  the  main  flow  was  shown  to  vary  directly  with 
the  turning.  With  increase  in  the  Mach  number  of  the  main- stream  flow, 
both  the  boundary-layer  thickness  and  the  deflection  increased. 


INTRODUCTION 

The  relative  lack  of  theories  explaining  the  behavior  of  boundary 
layer  when  a  lateral  curvature  of  the  main-stream  flow  exists  has  been 
especially  apparent  in  the  application  of  aerodynamics  to  the  design  of 
turbomachinery.  The  development  of  the  boundary  layer  in  such  cases  is 
strongly  influenced  by  the  corresponding  normal  pressure  gradient  toward 
the  center  of  curvature,  giving  rise  to  a  component  of  "secondary  flow" 
in  the  boundary  layer.  For  the  laminar  case,  most  of  the  published  work 
has  been  restricted  to  yawed  cylinders,  wings,  and  cones  (references  1 
to  4).  One  notable  exception  of  direct  application  to  the  design  of 
compressors  and  turbines  is  reference  5,  wherein  the  boundary  layer  on  a 
rotating  blade  is  analyzed.  For  the  turbulent  case,  a  general  but 
approximate  solution  of  the  momentum- integral  equations  based  on  an 
assumed  velocity  distribution  and  friction  law  is  obtained  in  refer¬ 
ence  6. 
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Aside  from  the  conventional  houndary-layer  approach  to  this  prob¬ 
lem,,  a  number  of  investigators  (references  7  to  9)  have  obtained  solu¬ 
tions  for  secondary  flow  arising  from  flows  of  varying  total  pressure 
or  varying  enthalpy  by  neglecting  the  influence  of  viscosity  but  admit¬ 
ting  the  existence  of  vorticity.  Although  such  procedure  obviously  does 
not  permit  satisfaction  of  all  the  boundary  conditions  (because  the 
order  of  the  general  differential  equations  for  the  flow  is  reduced) , 
the  results  so  obtained  give  a  fair  check  with  the  experimental  data, 
except  in  the  regions  close  to  the  wall.  Thus  the  indications  are  that, 
while  it  is  possible  to  obtain  a  fair  picture  of  three-dimensional  flows 
in  the  preceding  cases  by  neglecting  viscosity,  such  procedures  are 
inadmis sable  where  thin  boundary  layers  exist. 

The  object 'of  this  investigation,  which  was  conducted  at  the  MCA 
Lewis  laboratory,  is  to  solve  a  case  somewhat  analogous  to  that  of 
references  7  to  9  in  regard  to  main-flow  orientation  while  retaining  the 
usual  Prandtl  boundary-layer  assumptions,  and  thus  to  demonstrate  the 
influence  of  viscosity  and  compressibility  on  the  secondary  flow  in 
relatively  thin  boundary  layers. 

To  this  end,  the  main  flow  outside  the  boundary  layer  is  assumed  to 
follow  concentric  circular  streamlines  ih  planes  parallel  to  a  semi¬ 
infinite  flat  plate  and  to  be  uniform  in  a  direction  perpendicular  to 
the  plate.  The  edge  of  the  plate  is  placed  so  as  to  coincide  with  a 
radial  line  through  the  axis  of  the  main  flow  (fig.  l).  It  is  assumed, 
of  course,  that  in  the  established  main- stream  flow  the  effect  of  vis¬ 
cosity  is  neglibible.  In  the  simplification  of  the  equations  for  the 
flow,  it  was  found  convenient  to  follow  the  procedure  first  used  by 
Moore  (reference  4) . 

Although  the  solution  is  primarily  of  theoretical  interest,  it  may 
have  some  direct  bearing  on  the  understanding  of  secondary  flows  over  a 
portion  of  the  hub  between  the  compressor  guide  vanes  where  the  boundary 
layer  is  probably  still  thin  and  laminar,  and  where  the  main  flow  under¬ 
goes  relatively  little  turning,  often  along  a  circular  path.  Here  again 
the  analysis  would  probably  be  applicable  only  in  regions  where  boundary- 
layer  flows  from  the  blade  surfaces  and  other  "wall  effects”  were  not 
influential. 
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BASIC  EQUATIONS 

The  following  equations  describe  the  motion  of  a  compressible  vis¬ 
cous  fluid  in  vector  form  (reference  10 ): 

Conservation  of  mass 


2£+p(v.q)  =  0 


(1) 


Conservation  of  momentum 


p  m  -  Pf  "  VP  +  MS72*!  +  §  MV  (v*q)  +  2  £(vd)«v]  <1  +  (vn)x(7Xq)  - 


3  (v-q)(w) 


Conservation  of  energy 


P  §|  +  p(V-q)  =  A  +  V-  [j(kV)Tj 
where  the  dissipation  function  is  given  by 

A  =  d  |^2 V-  J(q-V)qJ  +  (VXq)2  -  2q-v(7-q)  -  |  (V-q)2 


Equation  of  state 


(2) 


(3) 


p  =  PRT  (4) 

(All  symbols  are  defined  in  the  appendix.) 

When  dealing  with  problems  involving  a  lateral  curvature  of  the 
potential  flow  it  is  convenient  to  use  an  orthogonal  coordinate  system 
illustrated  in  figure  2.  This  system  is  characterized  by  a  fixed  refer¬ 
ence  axis  x  of  arbitrary  curvature  c  =  c(x)  and  is  related  to  a 
Cartesian  system  as  follows: 


J 
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* 


*1  = 


cos  p  dx  +  z  sin  p 


x2  =  y 


)  (5) 


X3 


=  constant  +  z  cos 


sin  p  dx 


J 


where 


P  =  P(x) 

The  elements  of  length  at  (x,y, z)  in  the  direction  of  the  increas¬ 
ing  coordinates  are 

hj_dx  h2dy  h3dz 

and  their  values  may  he  obtained  by  use  of  equations  (5)  as 


% 


h^  =  (l  +  cz) 


h2  -  1 
h3  =  ! 

The  transformation  of  equations  (l),  (2) ,  and  (3)  from  the  Cartesian 
system  to  the  orthogonal  system  is  readily  obtainable  by  use  of  the 
expressions  for  the  elements  of  length  (see,  for  instance,  reference  11 
or  12) .  When  the  y-axis  is  assumed  perpendicular  to  the  surface  over 
which  the  flow  takes  place  and  the  usual  Prandtl  assumptions  about  the 
order  of  magnitude  of  various  terms  are  made,  the  equations  (l),  (2), 
and.  (3)  take  the  following  form  for  steady  flow  in  a  boundary  layer  where 
equations  (6),  (7),  and  (8)  correspond  to  equations  (l) ,  (2),  and  (3), 
respectively: 


2357 


2357 


NACA  IN  2658 


5 


* 

^  (Pu)  +  ^  [(l  +  cz)PvJ  +  ^  [(1  +  cz)pw  ►  =  0 

J 


(- 


du  du  Bu 


cz5x  +  v§y  +  w5^+l+czuw 


\  1  dp  d  /  du\ 

j  ~  ~  1  +  cz  5x  +  ciy  yP  §y J 


u  dw  3w  dw 


Tm'SL*  v  5^  +  w  §1  -  T 


Concomitant  with  (7a)  and  (7h)  is  the  result  that  ^  is  of  order  of 
magnitude  5  ~  0. 


u _  5E  SE  Se\  p  J  chi  5 


l  +  cz5x  +  v§y  +  *wr5z'/ - 1  + 


)  =  -  r+^ji?  +  |r  [0-  +  «0v]  + 


s  [d +  «h]  +  >*||?)  +(|f)  ]+  ^  (*  0)  (s) 

The  equation  of  state  (4)  remains  unchanged. 

Because  of  the  nature  of  the  problem  under  consideration,  several 
simplifications  in  the  preceding  equations  can  he  made.  First,  the  fixed 
reference  axis  x  will  he  chosen  to  coincide  with  a  streamline  of  the 
main-stream  flow  where  these  streamlines  are  assumed  to  he  arcs  of  con¬ 
centric  circles  in  planes  parallel  to  a  semi-infinite  flat  plate  (see 
fig.  l).  This  flat  plate  is  oriented  so  that  its  leading  edge  is  orthog¬ 
onal  to  the  reference  axis  xj  hence  the  leading  edge  can  be  defined  by 
x  =  0,  y  =  0.  Secondly,  if  the  velocity  of  the  fluid  in  the  main  stream 
is  designated  U]_,  it  follows  that  uj_  =  constant  along  any  given  stream¬ 
line;  that  is,  uj_  =  uq ( z ) . 

It  also  follows  from  the  nature  of  the  prescribed  main-stream  flow 
that  the  pressure  gradient  in  the  x-direction  and  outside  the  boundary 
layer  is  identically  zero.  In  view  of  the  order  of  magnitude  of  dp/dy, 
the  same  pressure  prevails  in  the  boundary  layer  as  in  the  main  flow, 
and  consequently  (where  the  subscript  indicates  partial  differentiation) 

Px  =  0  (9) 
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In  addition,  for  flow  with  circular  streamlines  it  may  he  shown 


Pz  ~  (  1  +  cz  )  plul^ 


The  boundary  conditions  on  equations  (6)  and  (7)  are 


u(x,«,z)  = 


w(x,»,z)  =  0 

u(x,0,z)  =  w(x,0,z)  =  v(x,0,z)  = 


where 


x  >  0 


and  at  the  leading  edge 


where 


u(0,y,z)  =  ux 
w(0,y,z)  =  v(0,y,z)  =  0 

y  >  0 


For  equation  (8),  usually  two  “boundary  conditions  on  temperature  are 
used  by  defining  its  value  at  y  =  0  and  y-*  co.  For  the  case  under 
consideration,  however,  it  will  be  necessary  to  define  temperature  only 
at  y  =  0. 


SOLUTION  OF  EQUATIONS 
Application  of  Vector  Potential  Concept 

In  order  to  satisfy  the  equation  of  continuity  (6)  and  to  eliminate 
one  dependent  variable  from  equations  (7),  the  concept  of  a  vector  poten¬ 
tial  is  employed  (see  reference  4).  To  this  end,  two  functions  ^  and 
cp  are  defined,  where 


2357 


MCA  TTT  2658 


7 


PU5+y 


(l  +  CZ)PW  =  Cpy 

(1  +  cz) pv  s  -(\|rx  +  cpz) 


J 


(12) 


It  may  be  verified  by  direct  substitution  of  equations  (12)  into  equa¬ 
tion  (6)  that  equation  (6)  is  identically  satisfied. 


Substitution  of  equations  (9),  (10),  and  (12)  into  equations  (7a) 
and  (7b)  yields 


t- 


-L 


+  9, 


y  \  p/x  Ty  \  p/z 

(l  +  cz) 


-  (tx  +  9Z)  \-2-L  + 


p/y  1  +  cz  Ty  \p 


9, 


*L 


p/y 


9y 

L  1 

ffz\ 

p(l  +  cz) 

+ 

X 

{l  +  cz) 

\p  / 

-  (tx+'Pz) 


9y 


c  9y  .  2 

- - rs  =  -  Cpn  U*i  2  + 

(1  +  cz)2  p  1  1 


p(l  +  cz) 

$ 


(13a) 


-  ?  ^y)2 


(13b) 


Incompressible  Case 

The  solution  for  the  incompressible  flow  (p  =  constant)  will  be 
presented  first.  For  this  case  it  is  necessary  to  consider  only  equa- 
tions  (13)  in  the  boundary  layer;  the  energy  equation  need  not  be  used. 

In  accordance  with  the  usual  procedure  in  the  solution  of  incom¬ 
pressible  boundary-layer  flows,  the  following  substitutions  for  the 
variables  are  employed: 


8 


MAC  A  TN  2658 


=  y  / 


ui 

Vx(l  +  zc) 


i|r  =  P  i/xpujlCi  +  zc)  F(t|) 

«p  =  pcx  (l  +  zc)  f^TTT)  G(tj)  ) 


(14) 


The  following  system  of  equations  is  obtained  by  substituting 
expressions  (14 )  along  with  their  appropriate  partial  derivatives  into 
equations  (l3)  (note  that  p/p  =  v) : 

(ex)2  [(14/c)  ar  * x]  - 

f,x)2  F"G  [|  +  %±-Hl  ^  -  F'  "  -  °  (15a) 


F'G'  -  + 


FG"  ‘  (l  -  F'2)  -  G' 1 • 


(l  +  zc) 

dull 

CU1 

dz  J 

(=x)2  G”G  [|  +  -  0 


(15b) 


It  is  interesting  to  note  that  setting 

u^  =  u^o(l  +  zc)n  (l6) 

in  (15a)  and  (l5b)  reduces  these  equations,  respectively,  to  the  forms 


(cx)2  G'F*  (n  +  1)  +  (ex)2  F"G  -  +  n  - 


2  -ci tin  l  3  „  \  _  i  t  _  o  (l7a) 


n(cx)2  (G*)2  -  (cx)2  G«g(|+  nj  +  F'G*  - 


FG" 


+  (l  -  F’2)  -  G* ‘ '  =0 

(17b) 
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Letting  n  =  -1  in  equation  (l6)  gives  the  expression  for  the  velocity- 
in  the  classical  irrotational  vortex  flow,  whereas  letting  n  =  1 
yields  the  expression  for  a  special  case  of  "wheel-type"  flow. 

Bow  F  and  G  and  their  derivatives  are  assumed  of  order  of  mag¬ 
nitude  of  unity ,  while  (cx)2«  1.  It  is  seen  that  cx  is  the  total 
turning  of  the  flow  outside  the  houndary  layer,  and  because  of  the  fore¬ 
going  restriction  it  must  remain  small.  It  is  further  assumed  that 

dup\ 

I  is  of  the  order  of  magnitude  of  unity  (which  should  be  the 

case  for  most  flows  of  interest  and  is  evidenced  in  part  by  equation  (16)). 
Beglecting  all  terms  of  the  order  of  (cx)2  in  equations  (15)  and  (17) 
reduces  these  equations  to  the  following  two  total  differential  equa¬ 
tions: 


(1  +  zc 
cux 


* 


■r 


FF"  +  2F' 1 »  =  0 


(18a) 


2F‘G»  -  FG"  -  2G*  "  +  2[l  -  (F‘)2J  =  0 

The  corresponding  boundary  conditions  are 

F=F'  =  G  =  G'=0atTj  =  0  ^ 

and  > 

F'  =  1  and  G*  =  0  as  ri  — »  » 

J 


(18b) 


(19) 


Equations  (l8)  and  (19)  show  that  for  small  turning  F  is  the 
usual  Blasius  function  and  that  equation  (l8b)  may  be  solved  for  G’ 

(or  G)  since  F  is  known.  Such  solution  was  obtained  by  relaxation 
procedure  inasmuch  as  F  and  its  derivatives  are  given  in  tabular  form 
only  (cf.  reference  13,  p.  121).  The  results  from  this  computation  are 
presented  in  figure  3  and  table  I. 


Boundary-layer  thickness.  -  If  the  boundary-layer  thickness  5  is 
defined  as  the  value  of  y  at  which  u  is  a  certain  percentage  of  u^ 
(usually  over  99  percent),  a  corresponding  value  for  T]  may  be  found 
from  the  solution  for  F  where  this  value  of  T)  is  designated  tj*. 
Consequently  the  expression  for  boundary-layer  thickness  may  be  obtained 
from  equations  (14)  as 


5  = 


(20) 
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or  nondimenslonally  in  terms  of  the  reference  length  x  as 

_1 

|  =  if  ^1  +  cz  Re  2  (21) 

On  a  given  streamline,  it  can  be  seen  from  equation  (20)  that  the 
boundary-layer  thickness  varies  directly -with  the  square  root  of  x. 

For  the  special  case  of  free  vortex  flow,  equation  (20)  becomes 

5  =  1*  (1  +  cz)  (22) 

and  it  follows  that  for  any  given  x  and  up^Q,  thickness  increases  in 
proportion  to  (l  +  cz)  or  that  the  boundary  layer  becomes  thinner  with 
an  increase  in  curvature  of  the  flow.  The  lines  of  constant  thickness 
are  obtained  by  considering  5  to  be  a  constant  in  equation  (22),  or 

_  constant 
(l  +  cz)^ 

Therefore,  the  curves  of  constant  8  are  spirals  of  lituus;  they  are 
shown  in  figure  4.  For  more  general  flows,  the  lines  of  constant 
boundary-layer  thickness  are  given  by 

u-^  (constant) 
x  _  (i  +  Zc) 

where,  of  course,  up  is  a  function  of  (l  +  zc). 

Flow  direction  at  plate  surface.  -  From  equations  (12)  and  (14)  the 
expressions  for  u  and  w  axe 


(23) 


u  =  Up  F1 
w  =  cxupG1 


Therefore,  if 


€  =  lim 
0 


w 

u 
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the  flow  direction  at  the  plate  is  given  "by 


a  =  arc  tan  e  (24) 

However,  since  G'(0)  =  F'(o)  =  0, 

€  =  CX  W)  (25) 

Equation  (25)  indicates  that 

T  _  £  T 
0,z  -  0,x 


inasmuch  as 


0) 


(y  -  0) 

The  ratio  G"(0)/F"(0)  was  evaluated  numerically  to  give  the  result 

€  =  -3.26  cx 

It  may  be  seen,  then,  that  at  the  plate  surface  the  flaw  is  deflec¬ 
ted  from  the  direction  of  the  streamline  toward  the  center  of  curvature 
of  the  main-stream  flow  and  that  this  deflection  is  proportional  to  the 
turning.  The  lines  of  constant  €  are  seen  from  equation  (25)  to  he 
those  of  constant  turning,  that  is,  radial  lines  (see  fig.  4). 


Effect  of  Compressibility 

Complexity  prevents  the  extension  of  the  previous  results  to  include 
the  effects  of  compressibility  in  the  same  manner  as  for  the  incompres¬ 
sible  case.  For  compressible  flow,  a  solution  can  be  obtained  only  in 
the  immediate  neighborhood  of  a  reference  streamline  on  which  the  Mach 
number  has  been  specified.  However,  it  is  possible  to  calculate  the 
-  Mach  numbers  on  various  chosen  streamlines  in  terms  of  the  Mach  number 

specified  on  an  initial  streamline.  Solutions  can  then  be  found  in  the 
neighborhood  of  the  chosen  streamlines  using  the  calculated  Mach  numbers. 
^  To  find  the  boundary-layer  flow  in  any  region  of  interest  is  therefore 

possible. 


12 


MCA  TN  2658 


Simplification  of  energy  equation,  -  Letting  0  =  CpT  and  assuming 
that  the  specific  heats  and  Prandtl  number  are  constant,  the  energy 
equation  can  he  written  in  terms  of  cp  and  as 


V* +  V*  -  (*= +  tpP‘V  -  rr^  “A +  +  «>>*  (fe) 


12 


{ 


p(l  +  cz)_ 


(26) 


For  Prandtl  number  equal  to  1  and  zero  heat  transfer  at  the  plate  sur¬ 
face,  equation  (26)  is  replaced  by 

°pt  + t  {(£/  -  [sirfsii] }-  T»°p  <27 

where  Tw  is  the  constant  temperature  of  the  plate. 

Reduction  of  equations  to  nondimensional  form.  -  The  following 
approximate  formula  (reference  14)  is  used  to  represent  the  relation 
between  the  viscosity  and  the  temperature: 


_J±_  =  A  _J_ 

^1,0  Tl,0 


(28) 


where  A  is  a  constant  suitably  adjusted  to  give  the  best  agreement 
with  the  actual  relation  over  the  temperature  range  considered. 

It  now  becomes  convenient  to  make  all  physical  quantities  in  the 
various  equations  nondimensional.  The  following  system  suggested  in 
reference  4  will  be  used;  all  quantities  on  the  left  are  to  be  con¬ 
sidered  relative  to  the  quantities  on  the  right,  which  are  measured 
outside  the  boundary  layer  and  along  the  reference  axis  x: 


u,v,w 

relative 

x,y,z,  i 
c 

relative 

p 

relative 

T 

relative 

to  u^o 

to  kl,oVpi,oui,o 

to  Pp^o 
to  u2i,o/2cp 
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p  relative  to  Pi,ou^l,0 

relative  to  (Ij^qA 

All  subsequent  equations  will  appear  in  nondimensional  form.  As  a 
consequence  of  tbe  use  of  the  nondimensional  form,  equation  (13a) 
becomes 


equation  (l3b)  becomes 


equation  (27)  becomes 


and  equation  (4)  will  be  written 

=  p 

pl,0  *1,0 

Transformation  of  Howarth.  -  The  following  transformation  of  vari- 
ables  indicated  by  Howarth  (reference  15)  is  now  utilized  to  make  the 
nondimensional  form  of  equations  (l3a)  and  (13b)  similar  to  the  incom¬ 
pressible  form  of  these  equations: 


(29) 
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pay 


"\ 
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«■ 


(30) 


As  a  consequence  of  transformations  (30),  equation  (29),  and  the 
fact  that  dp/dy  is  of  the  order  of  magnitude  5,  equations  (l3a), 
(l3b),  and  (27)  become,  respectively: 


¥y¥yx  +  $y¥yz  -  O^x  +  $z)¥yy  +  i  +  cz  ¥y$y 

-  %  {rhz  ^  *)+  (1  +  CZ)T™  (31s) 


^Y^YX  +  ®Y^YZ  "  (*X  +  ®zJ*YY  "  1  +  CS  ^Y^  +  CP  +  cZ^  Cp  U1 

*  +  W2  +  (its)  * 


Cfy)2] 

(31b) 

(32) 


The  same  procedure  as  for  the  incompressible  flow  case  is  now 
followed.  The  similarity  variable  and  corresponding  vector  potential 
components  are  defined  as 


v 


2357 


2357 


MCA  TN  2658 


15 


-\ 


u 


71  Y  yx(i  +  zc) 

¥  =  /xux(l  +  Zc)  F(t]) 


$  =  cX(l  +  Zc)  |/xu!(l  +  Zc)  G(tj) 


(33) 


By  use  of  these  definitions,  equations  (3l)  become 


(cX)2G'F» 


1  + 


^  S1]  ♦  <->2®"  [|  %2  - 1 


(l  +  Zc)  dul~l  FF" 


2cun  dZ 


-  F 


0 


(34a) 


FG 


F'G'  - 


-  +  ^1  +  -  F'2  ^ 


Q.M  + 


(cX)2G'2 


(l  +  Zc)  dul  r-1  »«  2 

l  cuj —  d z - 2“M 1  . 


2  '  2cux  dZ  2  1 


=  0 


(34b) 


Boundary  conditions  are  the  same  as  in  equations  (19)  with  the 
appropriate  change  in  variables;  that  is 


F  =  F'=G  =  G'=0  at  tj  =  0 


and 


F'  =  1  and  G'  =  0  as  T]-»oo 

_  Once  again  when  terms  of  the  order  of  magnitude  of  (cx)2  are 

eliminated  under  the  same  type  of  restrictions  as  in  the  incompressible 
case,  equation  (34a)  becomes  the  Blasius  equation  for  F.  However,  it 
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may  lie  seen  "that  equation  (34b)  cannot  be  solved  directly  as  a  total 
differential  equation  in  rf  since  Mi  is  a  function  of  Z.  In  fact, 
if  ui  is  specifically  prescribed  as  a  function  of  (l  +  cz)  (that  is, 
Ul  =  f(l  +  cz)  where  f  =  1  when  z  =  0) ,  then  Mp  may  be  expressed  as 


and  in  the  immediate  neighborhood  of  the  reference  axis  X,  Mp  «  M1?q 
and  may  be  considered  a  constant.  Incorporating  the  preceding  restric¬ 
tions  into  equation  (34b)  gives 

F'G'  -  G"'  -  l/2FG"+  (l  +  ^Ml,o)  f1  “  =  0  (36) 

This  equation  has  the  solution 

G*  =  (l  +  (37) 

where  G*  is  known  from  the  incompressible  case.  Curves  of  G*  (rj)  for 
various  values  of  are  plotted  in  figure  5. 

As  mentioned  previously,  the  flow  in  the  neighborhood  of  stream¬ 
lines  other  than  an  initial  reference  streamline  may  be  calculated  in 
the  same  manner  by  taking  the  reference  axis  along  the  new  streamline 
and  by  computing  the  appropriate  Mach  number  on  this  streamline  from 
equation  (35). 

Physical  coordinates.  -  Relations  (29),  (30),  and  (32)  may  be  used 
to  obtain  values  of  y  as  a  function  of  T)  for  a  fixed  value  of  x 


It  is  to  be  recalled,  of  course,  that  these  coordinates  are  the  non- 
dimensional  form  of  the  true  physical  coordinates. 

Interpretation  of  results.  -  From  equations  (18) ,  (36),  and  (38), 
the  effect  of  compressibility  may  be  seen  to  be  twofold,  namely,  changes 
occur  in  the  coordinate  y  as  a  function  of  r\  and  the  differential 
equation  for  G.  These  changes  have  a  corresponding  influence  on  the 
velocity  distribution,  the  boundary-layer  thickness,  and  the  flow  deflec¬ 
tion  at  the  plate  surface. 
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Since 


is  a  positive  function  of  T],  for  any 


given  value  of  rj  the  value  of  y  increases  vith  Mach  number,  as  shown 
by  equation  (38).  Consequently,  the  boundary-layer  thickness  5 
increases  with  Mach  number.  It  also  follows  for  a  chosen  Mach  number 


that 


I  -  Ml/2 

thereby  indicating  that  chi/By  increases  with  Mach  number  for 
7]  =  constant.  The  same  effect  is  known  to  exist  in  two-dimensional 
boundary-layer  flow. 

As  far  as  the  flow  deflection  at  the  plate  surface  is  concerned, 
the  changes  in  the  shape  of  G1  (see  fig.  5)  and  equation  (25)  indicate 
that  for  a  given  turning  the  deflection  toward  the  center  of  the  flow 
field  increases  with  Mach  number*  Consequently,  larger  secondary  flows 
in  the  boundary  layer  are  obtained  with  the  increase  in  Mach  number. 


1  +  lcrMi 


,o  i1  -  *'z) 


du 

dfj 


CONCLUSIONS 

The  following  conclusions  were  drawn  from  a  theoretical  investiga¬ 
tion  of  boundary-layer  flow  over  a  flat  plate  when  the  main- stream  flow 
outside  the  boundary  layer  follows  concentric  circular  streamlines  and 
the  total  turning  of  that  flow  is  small: 

1.  The  profile  of  the  velocity  component  in  the  direction  of  the 
flow  .outside  the  boundary  layer  is  identical  to  that  existing  in  two- 
dimensional  flow  over  a  flat  plate. 

2.  The  magnitude  of  the  velocity  component  in  the  direction  perpen¬ 
dicular  to  the  flow  outside  the  boundary  layer  and  in  a  plane  parallel 
to  the  plate  varies  directly  with  the  turning  of  the  main- stream  flow. 

3.  The  boundary-layer  thickness  varies  with  the  square  root  of  the 
distance  along  the  streamline  of  the  main- stream  flow  as  measured  from 
the  leading  edge  of  the  plate.  For  the  special  case  of  free  vortex 
flow,  the  lines  of  constant  boundary-layer  thickness  are  spirals  of 
lituus . 

4.  The  deflection  of  the  flow  at  the  plate  surface  from  the  direc¬ 
tion  of  the  main- stream  flow  is  directly  proportional  to  the  turning  of 
the  main-stream  flow.  The  lines  of  constant  deflection  are  therefore 
radial. 
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5.  Both  the  boundary-layer  thickness  and  the  deflection  at  the 
plate  surface  increase  with  Mach  number.  Consequently,  larger  secondary 
flows  in  the  boundary  layer  are  obtained  for  higher  Mach  numbers. 


Lewis  Flight  Propulsion  Laboratory 
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APPENDIX  -  SYMBOLS 


A 


cp 


The  following  symbols  are  used  in  this  report: 

constant  in  viscosity- temperature  relation 
specific  heat  at  constant  pressure  process 
curvature  of  x-axis,  dp/dx 


Dt 

V 

E 


G,G  J 

I 

h± 

k 

M 

P 

4 

B 

Re 

T 

U,v,w 

x,y,z 

a 


Eulerian  derivative 

vector  differential  operator 
energy 

functions  related  to  components  of  vector  potential 
inertial  field  force 

factor  used  in  transformation  of  coordinates  (i  =  1,  2,  3) 

thermal  conductivity 

Mach  number 

pressure 

velocity  vector 

gas  constant 

Reynolds  number  based  on  x,  (u^x/v) 
temperature 

velocity  components  in  curvilinear  coordinate  system 
coordinates  defined  by  equation  (30) 
curvilinear  coordinates 

boundary- layer  deflection  angle  measured  from  direction  of 
flow  at  plate  surface  to  direction  of  flow  outside 
boundary  layfer 
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P 

r 

8 

€ 


e 

A 


V 


p 

T0,x,  T0,z 


angle  between  Xj_-axis  and  tangent  to  x~axis 
ratio  of  specific  heats 
■boundary- layer  thickness 

measure  of  boundary-layer  deflection  near  plate  surface, 
tan  a 

similarity  variables 
value  of  T)  at  y  =  8 
enthalpy,  CpT 
dissipation  function 
absolute  viscosity 
kinematic  viscosity,  p/p 
density 

shear  stress  at  wall  in  x-  and  z-directions,  respectively 
components  of  vector  potential 


Cartesian  coordinates  (i  =  1,  2,  3) 


Subscripts: 

1  outside  the  boundary  layer 

0  on  reference  axis  x 

w  at  plate  surface 

All  other  subscripts  indicate  partial  differentiation* 
Superscript  primes  indicate  total  differentiation. 
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face  from  the  direction  of  a  circular  streamline  in  face  from  the  direction  of  a  circular  streamline  in 

the  main  flow  was  shown  to  vary  directly  with  the  the  main  flow  was  shown  to  vary  directly  with  the 

turning.  With  increase  in  the  Mach  number  of  the  turning.  With  increase  in  the  Mach  number  of  the 

main-stream  flow,  both  the  boundary -layer  thickness  main-stream  flow,  both  the  boundary -layer  thickness 

and  the  deflection  increased.  and  the  deflection  increased. 
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